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Waiting time jitter is a low-frequency phase jitter introduced whenever 
asynchronous digital signals are synchronized for multiplexing by pulse 
stuffing. It contains arbitrarily low frequency components and cannot be 
completely removed from the demultiplexed signals. 

In this paper the spectrum of waiting time jitter is derived, and sup- 
porting experimentally recorded waiting time jitter spectra are presented. 

A question of much engineering interest is at what rate waiting time 
jitter accumulates in a chain of multiplexer-demultiplexer pairs. A cal- 
culation based on the theoretical waiting time jitter spectrum shows that 
under reasonable conditions the rate at which the rms amplitude of the 
accumulated waiting time jitter grows is no greater than the square root 
of the number of multiplexer-demultiplexer pairs. Experimental data 
consistent with this upper bound are given. - 

I. INTRODUCTION 

When a number of lower speed digital signals are to be time division 
multiplexed to form a higher speed signal, it is necessary that the lower 
speed signals all be synchronous. If the sources of these signals are 
remote, this synchronization will be difficult to achieve. 

One way of achieving it is to distribute timing information from a 
master clock to the local clocks at the sources of the lower speed digital 
signals. In many situations, however, using a master clock and syn- 
chronizing the entire digital communication network is either not 
possible or not economically desirable. 

An alternative method of achieving synchronization is to use pulse 
stuffing synchronization. 1,2 The pulse stuffing synchronization scheme 
is reviewed in the next section. Basically, the idea is to bring the symbol 
rates of each of the incoming lines up to some common rate just before 
multiplexing by inserting dummy pulses. After demultiplexing, the 
dummy pulses are removed and each of the signals retimed to smooth 
out the gaps. 
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The retiming can never be done perfectly. There is always a phase 
jitter present on the demultiplexed signals. In Sections III and IV 
the spectrum of this degrading jitter is derived. Its accumulation in 
chains of multiplexer-demultiplexer pairs is studied in Section VI. 

IT. PULSE STUFFING SYNCHRONIZATION 

It is convenient to review pulse stuffing synchronization by describing 
the operation of a particular system. Since our experimental data 
were taken on the M12 multiplex, we have chosen this system to 
describe. 3 The reader should have no difficulty generalizing to other 
systems. The few details of the operation of the M12 that are essential 
for our discussion will be reviewed where necessary. 

The M12 is designed to multiplex four asynchronous DS-1 digital 
signals (nominally 1.544 Mb/s) into a DS-2 bit stream (nominally 
6.312 Mb/s). 4 For each of the four DS-1 digital data streams there is 
associated circuitry in the M12 designed to insert dummy pulses and 
bring its symbol rate up to a common one. 

A block diagram of one of these synchronizers is shown in Figure 1. 
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Fig. 1 — A block diagram of one of the four synchronizers in an M12. 



The data is written into the elastic store under the influence of the 
write clock, and read out under the influence of the read clock. The 
write clock is derived from the DS-1 signal. Its frequency (nominally 
1.544 MHz) will be denoted by / and its period, by / = l// n .* 
The clocking signal at the input to the inhibit gate is derived from 

* In making these definitions it is being tacitly assumed that there is no jitter 
on the write clock. This assumption will not always be reasonable since there may 
be jitter on the DS-1 signal from which it is derived. The effects of this jitter are 
studied in Section V. 
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a local clock. It is identical except for constant phase offsets with the 
corresponding clocking signals in the other three synchronizers. Thus, 
the data read out of the four synchronizers can be multiplexed by 
interleaving with no difficulty. 

The pulses on the clocking signal at the input to the inhibit gate 
are normally separated by four DS-2 time slots. Every twelfth pulse 
is separated from its successor by five DS-2 time slots to allow framing 
and control pulses to be inserted onto the DS-2 signal. The average 
frequency of this clocking signal will be denoted by / x and its average 
period by t x = 1/fi . Nominally, 

/, = (12/49) -6.312 = 1.5458 MHz. 

For the remainder of this paper, the jitter on this clocking signal will 
be neglected, and it will be approximated as an unjittered clocking 
signal with frequency /j .* 

By design the frequency /, is greater than /„ . If all the pulses at the 
input to the inhibit gate were allowed to pass through, the read clock 
would eventually overtake the write clock. Since it is necessary to 
write before reading, some pulses must be inhibited. 

The inhibiting of pulses is done by the phase comparator and as- 
sociated control circuitry. The output of the phase comparator, <t>s(t), 
is the phase difference between the read and write clocks (or, equiv- 
alent^, the jitter on the output data stream of the synchronizer). 
The control circuitry monitors 4> s (t) and as soon as possible after <f> s (t) 
crosses some predetermined threshold A inhibits a clocking pulse, or, 
as it is more commonly phrased, stuffs a (data) pulse. 

A typical graph of <f> s (t) is shown in Figure 2. By convention <f> s (t) 
is defined as positive if the read clock precedes the write clock. This 
seemingly backwards definition gives <f> s (t) positive slope and negative 
jumps as seems to be traditional. 5 ,e Since the write clock must always 
precede the read clock, the threshold A must be set so that <f>s(t) is 
always negative. 

The positive slope of 4> s (t) at nonjump times reflects the fact that 
in the absence of inhibit pulses the read clock is gradually overtaking 
the write clock. Whenever a pulse is stuffed, if>s(t) drops by one slot 
(by definition one slot equals ti seconds). Stuffing does not occur as 
soon as <t>s(t) crosses the threshold. It is necessary to wait until a stuffing 



* The techniques of Section V could be used to remove this assumption also. It 
is not particularly interesting to do so, however, since little error results from making 
it. The jitter being neglected is systematic and of relatively high frequency. It is 
easily removed at the demultiplexer. 
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Fig. 2 — A typical graph of the output of the phase comparator. 



opportunity occurs. These opportunities, which are fixed by M12 
format considerations, occur at regular intervals 1176 DS-2 time slots 
apart. Let t, n denote this interval. Then, 

t m = 1176 -(12/49) ty 

= 288 U . 

The maximum possible stuffing rate is / m = l/t m . In the M12 

/ m = 6.312/1176 MHz 

= 5.367 kHz. 

The actual average rate at which pulses are being stuffed is /. = /i — /o ■ 
The stuffing ratio 

P = 1./U , 

which must be between and 1, is a critical parameter. Nominally, 
in the M12 

p = LIU = (U - U)/U 

_ (12/49) -6.312 - 1.544 
6.312/1176 

= 0.3346, 

or slightly more than 1/3. 

It will be convenient for the remainder of this paper to measure 
time in stuffing opportunities and frequency in cycles per stuffing 
opportunity. All numerical values for times and frequencies given 
from here on are to be assumed to have these units unless other units 
are explicitly stated. Conversion to seconds and Hertz is easy through 
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the relations 



and 



1 stuffing opportunity = 0.186 ms 
1 cycle per stuffing opportunity = 5.367 kHz. 



It has been traditional to consider <j> s (t) as being the sum of two 
jitter waveforms — the first being stuffing jitter and the second waiting 
time jitter. Stuffing jitter is defined as the jitter that would be present 
if stuffing could occur on demand, and waiting time jitter is usually 
defined as the low-frequency jitter present in 4> s (l) because in actuality 
there is a "waiting time" between stuff demand and stuff execution. 
Since these definitions are hard to make precise, we shall not attempt 
to make any such distinctions and shall merely call <f> s (l) the waiting 
time jitter waveform. 

After each of the four DS-1 signals is synchronized, it is multiplexed 
with the other three, transmitted over an appropriate facility and 
demultiplexed by an M12. Stuffed pulses on each of the lines are removed 
in an associated desynchronizer. 

A block diagram of a desynchronizer is shown in Figure 3. Information 
telling the demultiplexer where stuffs have been made is carried on the 
DS-2 signal along with the data. Thus, it is possible for the write clock 
in the desynchonizer to exactly duplicate the read clock in the syn- 
chronizer. It has gaps where stuffs were made. 

The read clock on the elastic store in the desynchronizer is obtained 
from a voltage controlled oscillator (VCO). The output of the phase 
comparator, which is proportional to the phase difference between 
the read and write clocks, is filtered and fed back to the VCO in such 
a way as to keep the read clock from overtaking or falling too far 
behind the write clock. 
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Fig. 3 — A block diagram of one of the four desynchronizers in an M12. 
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The desynchronizer acts as a phase-locked loop and smooths the gaps 
created by stuffing pulses. Ideally, the output of the desynchronizer 
would be identical with the input to the synchronizer. Unfortunately, 
all phase-locked loops have a low-pass nature and some low-frequency 
jitter degrades the output. 

Let 4> D (t) denote this jitter. Then 

<t>n(t) = h(t)*<t> s (t), 

where h(i) is the overall transfer function of the phase-locked loop in the 
desynchronizer and * denotes convolution. In the M12 h(t) falls off 
at 12 dB per octave after cutoff at about 0.12 cycle per stuffing op- 
portunity (644 Hz).* 

Often <f> s (t) and <f> D (t) are both called waiting time jitter (dis- 
tinguishable only by context). To avoid any confusion, we shall only 
call (j>s (0 waiting time jitter and shall always call <f> D (t) filtered waiting 
time jitter. 

Because of the low-pass nature of h(t), the low-frequency components 
of waiting time jitter are present in filtered waiting time jitter and 
degrade the demultiplexed signal. It is not difficult to demonstrate 
that it is quite possible for 4> s (t) to contain significant low-frequency 
power. In Figure 4a, 4>s(t) is drawn assuming p = 1/2. The lowest 
frequency component in this waveform is at the relatively high frequency 
of 1/2 cycle per stuffing opportunity. But, it is unrealistic to assume 
p = 1/2 exactly. Figures 4b and 4c show <f> s (t) for p slightly greater 
than 1/2 and p slightly less than 1/2. In both cases there is a strong 
low-frequency envelope. In general, whenever p is close to a simple 
rational number (a rational number with a small denominator), but 
not exactly equal to it, <j> s (t) will have appreciable low-frequency 
power. 

This idea is developed in papers by Kozuka 5 and Matsuura, Kozuka, 
and Yuki . Formulas are given for the peak-to-peak amplitude of the 
low-frequency envelope present on the waveform <j> s (t) when p is close 
to, but not exactly equal to, a simple rational number. Questions 
arise, however, as to how close is close enough and how the trade off 
between approximation accuracy and simpler rational numbers is to 
be made (for example, should 0.45 + be considered close to 1/2 or to 
9/20?). 

An alternate approach is to calculate the power spectrum of waiting 
time jitter as we do in the next section. 



* The Ml 2 has been recently redesigned and now has a lower cutoff frequency. 
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Fig. 4 — In the upper portion of this figure p equals 1/2 exactly, and there are 
no low-frequency components in <P S (I)- In the other portions of this figure p is close 
but no) exactly equal to 1/2, and there are strong low-frequency components. 



III. THE SPECTRUM OF WAITING TIME JITTER 

3.1 Theoretical 

In Appendix A the power spectrum S s (f) of 0. S (O is calculated. 
Briefly, the procedure used is to (/) find an equation describing waiting 
time jitter waveforms, (ii) introduce initial condition random variables 
into this equation in such a way that a stationary ensemble of waiting 
time jitter waveforms is defined, (Hi) compute the covariance of the 
waiting time jitter random process, and (iv) Fourier transform this 
covariance to obtain the power spectrum S s (j). 

It is possible to find the power spectrum of waiting time jitter without 
using random process theory. A procedure similar to one used by 
Iwereon 7 to find the spectrum of quantizing noise in a delta modulator 
with unequal positive and negative step sizes could be used. (The 
similarity of the underlying mathematics in these two at first seemingly 
unrelated problems is remarkable.) Our method of finding the spectrum 
(that is, via random processes) does not require any approximations. 
The method used by Iwerson requires approximating pulses by delta 
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functions. This approximation is unimportant for his problem, but of 
more consequence in the waiting time jitter problem. A further 
advantage of proceeding via random process theory is that the arguments 
needed are somewhat simpler and more concise, especially those that 
are given in Appendix B of this paper (corresponding with those in 
Appendix D of the paper by Iwerson). 

The result of Appendix A is (remember, frequency is being measured 
in cycles per stuffing opportunity) 

S s (f) = sine 2 j-QiJ) + g GfJ'W/ - n ) + W + n ))> W 

where 

, sin irf 
sine / = — t~ , 

TV) 

Q(J) = Z (as)Vq> S( J ~ pn) + rep 5(/ + pn)) - (2) 

repX(/)= £ X(j - k) for any function X(j), 

k = -oo 

and <5(-) is the Dirac delta function. It is convenient to define 

S s , A (j) = sinc 2 /-Q(/) (3) 

and 

Ss.bU) = E Gfjw - ») + W - n »- < 4 > 

With this notation 

S s (1) = S S , A (1) + SsAf)- 

Graphs of sine 2 /, Q(f), and S 8 , B (f) are given in Figure 5. Only the 
lines associated with the first, second, and third terms of the defining 
sum for Q(f) are shown in the graph of Q(f). The lines labeled 1" are 
lines introduced by 

y^repaff-p-l). 

Lines labeled 1 + are lines introduced bv 



1 

2tt-1 



repfitf + p-1). 
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Fig. 5 — The components of the spectrum of waiting time jitter. 



In general, n -lines and w + -lines are lines introduced by 

1 



and 



2m 



2-jrn 



rep 8(1 — pn) 



rep 8(1 + pn) 



respectively. This notational scheme is similar to the one used by 
Iwerson. 

The interesting information about the spectrum S s (l) is carried 
by Q(D- S s ,u(l) contains only relatively high-frequency components. 
The function sine 2 / is just an envelope. The low-frequency content 
of S s (1) is determined by p through Q(l). 

If p is irrational, none of the spectral lines will coincide. When p 
is rational, the spectral lines do eventually coincide and it is possible 
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to replace the infinite sum in the denning equation for Q(f) by a finite 
sum. Iwerson shows that if p = p/q where p and q are relatively prime, 



Q(f) = ^ § esc 2 (%)rep *(/ - **) + ^ rep *(/). 



(5) 



3.2 Experimental 

The spectra of waiting time jitter waveforms produced by an M12 
were recorded for a number of stuffing ratios. Three representative 
spectra are shown in Figures 6, 7, and 8. The stuffing ratios in these 
three figures are respectively 0.186, 0.372, and 0.333. The stuffing ratio 
was varied by varying the DS-1 symbol rate. 

The vertical dB scale on the graphs is only approximate. In addition, 
no attempt was made to have dB correspond to any particular amount 
of jitter power, and therefore, readings taken from it are only relative. 

In Figures 6 and 7 the stuffing ratio is not close to any simple rational 
number and the spectral lines are all distinct. To record Figure 8, 
the stuffing ratio was set as close as possible to 1/3. The expected 
coincidence of lines is quite apparent. 

The weak line at about 0.08 cycle per stuffing opportunity in Figure 8 
is not present in the theoretical spectrum. Its origin can be traced to 
an improper (only in the sense that it deviates from the usual model) 
functioning of the synchronizer in the M12. This improper functioning, 
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Fig. 6 — An experimental waiting time jitter spectrum with p = 0.186. 
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Fig. 7 — An experimental waiting time jitter spectrum with p = 0.372. 

termed multiple stuffing, is also apparent when looking at oscilloscope 
traces of waiting time jitter waveforms. Occasionally stuffs will be 
made when they quite obviously should not be (for instance, twice 
in a row with p < 1/2). The actual Ml 2 used to record the spectra 
of Figures 6, 7, and 8 was an experimental model having an asynchronous 
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Fig. 8 — An experimental waiting time jitter spectrum with p = 0.333. 
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delay line 8 elastic store rather than the usual commutator type (see 
page 617 of Ref. 4). It was used because it exhibited slightly less multiple 
stuffing. 

In Tables I, II, and III the powers of the lines in the experimental 
and theoretical waiting time jitter spectra are tabulated for the stuffing 
ratios 0.186, 0.372, and 0.333 respectively. To obtain the theoretical 
powers listed in Tables I and II, equation (2) for Q(/) was used. The 
rational stuffing ratio formula (5) for Q(j) was used in computing 
the theoretical powers listed in Table III. 

The entries in the columns headed uncorrected experimental power 
were simply read off Figures 6, 7, and 8. For ease of comparison cor- 
rected experimental powers are also listed. These entries were obtained 
by multiplying the corresponding uncorrected powers by 0.817 (to 
correct the vertical scaling) and subtracting 42.0 (to add an absolute 
reference). The constants 0.817 and 42.0 were chosen by a least squares 
procedure. 

The differences between the theoretical and corrected experimental 
powers are listed in the last column of Tables I, II, and III. The agree- 
ment is to within 0.5 dB except in a few cases where discrepancies are 
easily explained in terms of the low frequency and small amplitude 
limitations of the spectrum analyzer. 

IV. THE RMS AMPLITUDE OF FILTERED WAITING TIME JITTER 



4.1 Theoretical 
Since 



**(*) = h(t)*<j> s (t), 



Table I — Spectral Line Powers for p = 0.186 





Normalized 


Power (dB) 






Experimental 


Experimental 




Line 


Frequency 


Theoretical 


Uncorrected 


Corrected 


Difference 


1- 


0.186 


-16.5 


31.0 


-16.6 


-0.1 


2" 


0.372 


-24.1 


22.3 


-23.7 


0.4 


3~ 


0.558 


-30.5 


14.1 


-30.4 


0.1 


4" 


0.744 


-38.2 


5.6 


-37.4 


0.8 


5- 


0.930 


-52.5 


— 


— 


— 


1 + 


0.814 


-29.3 


15.2 


-29.5 


-0.2 


2+ 


0.628 


-28.6 


16.1 


-28.8 


-0.2 


3+ 


0.442 


-28.5 


16.3 


-28.6 


-0.1 


4+ 


0.256 


-29.0 


15.8 


-29.0 


0.0 


5 + 


0.070 


-30.0 


14.3 


-30.3 


-0.3 


Ss.D 


1.000 


-30.6 


13.5 


-30.9 


-0.3 
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Table II — Spectral Line Powers for p = 0.372 





Normalized 


Power (dB) 






Experimental 


Experimental 




Line 


Frequency 


Theoretical 


Uncorrected 


Corrected 


Difference 


1- 


0.372 


-18.0 


29.4 


-17.9 


0.1 


2- 


0.744 


-32.2 


11.8 


-32.3 


-0.1 


3- 


0.116 


-25.7 


20.1 


-25.5 


0.2 


4~ 


0.488 


-31.8 


11.8 


-32.3 


-0.5 


5~ 


0.860 


-46.0 


1.2 


-41.0 


5.0 


1 + 


0.628 


-22.6 


24.0 


-22.4 


0.2 


2+ 


0.256 


-22.9 


23.8 


-22.5 


0.4 


3 + 


0.884 


-43.3 


2.0 


-40.3 


3.0 


4+ 


0.512 


-32.1 


11.6 


-32.5 


-0.4 


5+ 


0.140 


-30.2 


14.0 


-30.5 


-0.3 


03, B 


1.000 


-24.6 


21.3 


-24.6 


0.0 



the spectrum S n (f) of filtered waiting time jitter is given by 
S D (j) = \H(j)\ 2 S 8 (f) 

where H(f) is the Fourier transform of h(t), that is, 
H(j) - 5\h{t)) 



-/: 



h(t)e 



-i2T ft 



dt. 



Only the low frequency components of S s (f) will not be significantly 
attenuated by | H (/) | 2 . 

The power P»(p) in <p n (t) (that is, the rms amplitude squared of 
<M0) is given by 

pm = r s D (f) df 

J -co 

= r 1 inn 1 2 s s a) df 

J -co 

= J" I H(j) I 2 S s . A (f) df + p I H(J) | 2 S s .„(f) df 

4 P n .M + P D . B (p). 

Typically, h(() cuts off at frequencies much smaller than 1 cycle per 
stuffing opportunity where the lowest frequency spectral line in Ss.nif) 
is located, and, therefore, P, KB (p) is quite small in comparison to 
Pd,a(p)- Since Pd, a (p) is symmetric about 1/2 (that is, P d ,a(p) = 
Pd,a(1 ~ p))> >t is to be expected that P d (p) will be essentially sym- 
metric also. 
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Fig. 9 — A theoretical graph of the rms amplitude of filtered waiting time jitter as 
a function of p with H(f) assumed to have a double pole at 0.12 cycle per stuffing 
opportunity. 



A computer program was written to plot P d (p) for various filter 
transfer functions H(f). Two sample graphs appear in Figures 9 and 10. 
The assumed transfer functions were respectively 



and 



H(j) = 



H(f) = 



0.12 



if + 0.12 



0.06 



\jf + 0.06, 

The first transfer function is a good approximation to that of the 
filter in the M12*. The cutoff frequency of the second transfer function 
is 1/2 that of the first. In both graphs the predicted symmetry is quite 
apparent as are the expected peaks at rational stuffing ratios. The 
peaking is more pronounced with the sharper filter. 

For both filters there appears to be significant benefit to be gained 
by choosing p intelligently. Actually, as will be discussed more fully 
in Subsection 5.2, this conclusion is misleading because the effect of 
jitter on the input DS-1 signal has not been considered. 



* As noted in a previous footnote, the phase-locked loop in the M12 has been re- 
cently redesigned and now has a lower cutoff. 
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Fig. 10 — A theoretical graph of the rms amplitude of filtered waiting time jitter 
as a function of p with H (/) assumed to have a double pole at 0.06 cycle per stuffing 
opportunity. 

4.2 Experimental 

An experimental recording of Pi>(p) versus p for p t [0, 1] appears 
in Figure 11. As in the experimental spectra of Figures 6, 7, and 8 
the dB scale is only approximate and relative. With this fact taken 
into account, the agreement with Figure 9 is good. 

The vertical lines at the tops of the peaks are caused by the low- 
frequency cutoff of the true rms voltmeter used in making Figure 11. 
The unit used cutoff at 10 -4 cycle per stuffing opportunity (0.5 Hz). 



V. INPUT JITTER 

5.1 Generalized Spectrum 

It is possible to generalize the spectrum derivation of Appendix A 
to allow for jitter already present on the signal at the input to the 
synchronizer. This jitter will indeed usually be present. If the signal 
at the input to the synchronizer is coming from a moderately sized 
Tl line, it will be degraded by significant repeater jitter.""" In addition, 
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it will be degraded by filtered waiting time jitter if it has already passed 
through a syndes (synchronizer-desynchronizer pair). 

Let 0,(0 denote the jitter on the signal at the input to the synchronizer 
and <j) S (t) and <t>/,(t) again denote the jitters on the outputs of the syn- 
chronizer and desynchronizer. All of these jitters will be assumed to 
be defined with respect to an unjittered reference signal and to be 
positive if their associated signals are ahead of it. In Section III the 
output of the phase comparator and the jitter on the output of the 
synchronizer were identical, and we denoted them both by lS (O- When 
there is input jitter, these two jitters are no longer identical and a new 
symbol must be introduced for one of them. We have chosen to still 
call the jitter on the output of the synchronizer 0s (0- The symbol 
<t>spc(t) will be used for the output of the synchronizer's phase com- 
parator. We have 

+s*c(0 - 0.s(O - 0/ CO- 
In Appendix B it is shown that with input jitter accounted for the 
spectrum of waiting time jitter is given by 

8,® - sine 2 /•(?(/) + £ (^J'W - ») + «(/ + n)) 



= S a . A (f) + S„. B (j) 



+ sine 2 /-rep £,(/) 

+ Ss.cd), (6) 




0.4 0.5 0.6 
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Fig. 11 — An experimental graph of the rms amplitude of filtered waiting time 
jitter as a function of p. 
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where 



Q(1) = £ (^JV e P Z »U ~ n "> + re P Z ^ + np)) 

ZJS) - ffWO) 

z n (t)e- i2 *" dt, 



(7) 



=/: 



«„(*) = #(exp {-jam(0/(O - 0/(0))}}, 

£,(/) is the spectrum of the input jitter, and E denotes expectation. 

The input jitter fait) will usually only have significant power in 

frequencies that are quite small in comparison to 1 cycle per stuffing 

opportunity. When this is the case, the approximation (see Figure 12) 





■ro^. 
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Fig. 12 — Approximating sine 2 / rep S,(f) by S,(f). 
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sine 2 /-rep S t (f) = sine 2 f-S t (f) 
= Sj(j) 

can be used to simplify (6) to 

Sa(f) = sine 2 /■£(/) + £ (^W ~ ») + *(/ + »)) + AW- (8) 

The input jitter makes its presence felt in two ways. The first way 
is simply in the appearance of the feedthrough term Si(j) (or, more 
precisely, sine 2 /-rep &/(/)). The second way is in the smearing of 
spectral lines. The function Q({) has the same functional form as Q(j), 
but the envelope functions Z„(/) replace impulses. 

The envelope functions Z H (f) are difficult to compute in general. 
Evaluating them requires evaluating the same type expressions that 
must be evaluated to find the spectrum of exponentially modulated 
carriers. This latter problem has been a traditionally difficult one. 
Exact evaluations for the special cases of Gaussian input jitter and 
sinusoidal input jitter and approximate evaluations for small amplitude 
input jitter are possible. 

5.1.1 Gaussian Input Jitter 

If tf>i(l) is Gaussian with zero mean and covariance Ci(t), then the 
random variable 0/(0 — 0/(0) is Gaussian with zero mean and co- 
variance 2(C/(0) — C/(/)). By the well known formula for the char- 
acteristic function of a Gaussian random variable, 

z n (t) = E{exp \-j2 v n(4>i{t) - 0/(0))}} 

= exp{-(2,m) 2 (C/(0) - CM)- 

Using the infinite series expansion for the exponential function, we have 

z n {t) = exp j-(27m) 2 C/(0)} exp \{2^n) 2 C,{t)} 

= exp {-(27rn) 2 C/(0)} £®?^Cj(0. 

*=o Kl 

Therefore, 

= exp {-(27m) 2 C/(0)} ± &%f- ff { C)(t) } . 

= exp{-(27m) 2 C/(0)} 8(f) 
+ exp | - (27rn) 2 C/(0) } £ ^jf- S,(fl* ■ ■ ■ *S,(f) . (9) 

t_l ft I • v ' 

k terms 
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The envelope Z n {j) is composed of an impulse at the origin (the 
carrier line in exponential modulation theory) and an infinite sum of 
Mh order convolutions of the input jitter spectrum with itself. 

5.1.2 Sinusoidal Input Jitter 

If the input jitter is sinusoidal, exact results are also possible. Sinus- 
oidal input jitter can be modeled by the random process 

<f>j(t) = /3sin (2vFt + 0), 

where and F are constants and is a random variable distributed 
uniformly on [0, 2ir). [Introducing is necessary to randomize the 
epoch and make0 7 (-) stationary.] 
By direct substitution, 

z n (t) = E{ex V \-j2MMt) ~ 0/(0))}} 

= E\exp \-j2m0sm (2-n-Ft + 0)} exp {j27m/3sin 0}} 
= #{exp {j'27rn/3sin {-2irFt - 0)} exp \j2irnp sin 0}}. 
Using the expansion 

e ,xainx = £ J k ( X )e ikI , 

k=-oo 

where the J k (\) are Bessel functions of the first kind (see, for example, 
equation 7.53 of Ref. 12) in this equation, we have 



00 



z n (t) = B< S J k (2Tm$)e m2TFt + e) • E Ji(2*npW 

U=-oo i--oo 

= £ J k (2T r np)J l (2Tnf3y ik2 * F, E{e'"- k) °} 

fc.f— so 

= £ J 2 k (2Tm®e- ik2 "" . 
Therefore, 

Z n (f) = Z Jl&*n(3) W + W- (10) 

* = -oo 

The transform Z n (f) consists of a string of impulses stretching from 
— oo to +oo and separated from each other by F cycles per stuffing 
opportunity. 

5.1.3 Sin all Sign at Ap proxim ation s 

Small signal approximations for the Z„(f) can be derived easily 
using the power series expansion for an exponential function. We have 
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z n (t) = £{exp {-jftrn(<fc(0 - 0/(0))}} 
= 1 - j2imE{<f, I (t) -0/(0)} 

+ ^ (-i2,m) 2 ^{(0 J (O - 0/(O)) 2 } + • • • 

= 1 - KWtf {(0/(0 - 0/(O)) 2 } + • ■ • 
= 1 - (2^n) 2 (C/(0) - C,(0) + • • • , 



and 



Z n (j) = fffc(0] 

= «(/) - (2*n)\C,(Q) 5(/) - &(/)) + ■ • ■ . (11) 

If n and the input jitter are small enough, it will be reasonable to 
approximate Z n (f) by the first few terms of this expansion. 
The simplest such approximation is 

Zn(J) = «(/)• (12) 

Taking one more term gives the better approximation 

Z n (f) = *(/) - (27m) 2 (C 7 (0) 8(f) - S t (f)) 

= (1 - (2,m) 2 C/(0)) 8(f) + (27m) 2 >S/(/). (13) 

Both of these approximations have the desirable property of preserving 
total power. 

The power in the term (2ini) 2 Sj(f) can be used as a guide for deter- 
mining which of these approximations (if either) is reasonable. For 
all n, the total power in Z„(/) is 1. If 

(27m) 2 P Sj(j) df 4 (27m) 2 P 7 

J— to 

is very small in comparison to one (say < 0.01), the approximation 

ZJj) = 5(f) 

will be adequate. For moderate power (say 0.01 < (27r??) 2 P, < 0.1), 
the second approximation will still be reasonable. 

For a given input jitter power, either of these approximations is 
better for small n than large. For weak input jitter 

Z n (f) = 8(f) 

will be reasonable for small ??, but not for large. Thus, the effect of 
weak input jitter is to smear the high order rc-lines. [The feedthrough 
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term will be negligible because if Z n (j) = 8(f) is reasonable for even 
small n, sine 2 /-rep £,(/) will be negligible in comparison to sine 2 /•$(/)]. 

5.2 Filtered Waiting Time Jitter Power with Gaussian Input Jitter 
When the effect of input jitter is considered, the power in filtered 
waiting time jitter is again given by 



Up) = f I H(j) | 2 s a (fi df, 

J — to 



but with Ssd) as defined by equation (6). Let P d .a(p), Pd.b(p), and 
P DiC denote the components of P d (p) due to S s .a(1), S s , B (f), and 
S a .c(j) respectively. The power P DiC is a constant independent of p. If 

sine 2 /-rep S x (j) = £,(/), 

then 

p D .c= r \Hd) i 2 s^df, 

J —oo 

which is the power of the jitter that would be present at the output 
of the syndes (synchronizer-desynchronizer pair) if its only effect 
were to filter the input jitter. The jitter power added to the bit stream 
by pulse stuffing is to within this approximation 

APM = P D .M + P d ,b(p)- 

Computer-drawn graphs of AP D (p) for Gaussian input jitter appear 
in Figures 13 through 16. The transfer function H{j) was assumed to 
equal that of Figure 9 for Figures 13 and 15 and that of Figure 10 for 
Figures 14 and 16. For Figures 13 and 14 

Sl(j) = 2^ liO/ (0.1)' + f ' 
and for Figures 15 and 16 

Sl(j) = %r W (0.1) 2 '+f 

Both of these spectra are RC with a cutoff at 0.1 cycle per stuffing 
opportunity. The rms amplitude of the jitter of the first is 1/10 slot, 
and of the second, 1/4 slot. The first jitter spectrum is felt to be repre- 
sentative of what will typically be encountered in the field. 

A comparison of Figures 9, 13, and 15 and of Figures 10, 14, and 16 
shows that the effect of input jitter is to erode the peaks and valleys 
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Fig. 13 — A theoretical graph of the jitter power added by pulse stuffing as a 
function of p. The jitter on the input to the syndes is assumed to be Gaussian with 
an rms amplitude of 1/10 slot and an RC spectrum with a corner frequency of 0.1 
cycle per stuffing opportunity. The transfer function H(f) is assumed to have a 
double pole at 0.12 cycle per stuffing opportunity. 
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Fig. 14 — A theoretical graph of the jitter power added by pulse stuffing as a 
function of p. The jitter on the input to the syndes is assumed to be Gaussian with 
an rms amplitude of 1/10 slot and an RC spectrum with a corner frequency of 0.1 
cycle per stuffing opportunity. The transfer function H(f) is assumed to have a 
double pole at 0.06 cycle per stuffing opportunity. 
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Fig. 15 — A theoretical graph of the jitter power added by pulse stuffing as a 
function of p. The jitter on the input to the syndes is assumed to be Gaussian with 
an rms amplitude of 1/4 slot and an RC spectrum with a corner frequency of 0.1 
cycle per stuffing opportunity. The transfer function H(f) is assumed to have a 
double pole at 0.12 cycle per stuffing opportunity. 

of the no input jitter graphs. It is felt that this erosion will occur for 
almost all input jitter and not just for Gaussian input jitter with a 
spectrum of the above form. Thus, the choice of p is not so critical 
as it appears to be in Figures 9 and 10 when, as is usually the case, 
input jitter is not negligible. Notice, however, that even with non- 
negligible input jitter, there is still an advantage in using a narrow 
phase-locked loop. The graphs of Figures 14 and 16 are typically 3 dB 
below those of Figures 13 and 15. 

5.3 A Bound on the Power in Filtered Waiting Time Jitter 

By neglecting the effect of the filter in the desynchronizer, it is 
possible to obtain a simple bound on the power in filtered waiting time 
jitter. This bound is quite weak, but will be useful in the next section 
to bound the accumulation rate of filtered waiting time jitter in chains 
of syndes. 

Let 

H max = max | H(j) I. 



Then 



P„ ^ H* m „ P s 
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where 

Ps = P S a (J) dj 

J -oo 

= r ssM dj + r s a . fl a) d/ + r s a . C 0) # 

J-M »"-00 "-00 

= ■* S.A I *S,A T *S.C • 

From equations (32), (31), and (33) of Appendix B 

Ps,x - 1/12, 

P S . B = P 2 /12 £ 1/12, 
and 

Ps.c = "/ t 
where P 7 is the input jitter power. Therefore 

P D < HLx(l/6 + P;). 

If there is no peaking, 

tf m!lx = 7/(0) = 1 
and 

Pn ^ 1/6 + Pi . 



(14) 



05) 
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Fig. 16 — A theoretical graph of the jitter power added by pulse stuffing as a 
function of p. The jitter on the input to the syndes is assumed to be Gaussian with 
an rms amplitude of 1/4 slot and an RC spectrum with a corner frequency of 0.1 
cycle per stuffing opportunity. The transfer function H(f) is assumed to have a 
double pole at 0.06 cycle per stuffing opportunity. 
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VI. WAITING TIME JITTER ACCUMULATION 

6.1 Theoretical Bound 

A question of much engineering interest is how fast filtered waiting 
time jitter accumulates in tandem connections of syndes. Although 
the bound (15) is quite weak, it does guarantee that when there is no 
peaking in any of the transfer functions of the desynchronizers in the 
chain, the rms amplit ude o f the jitter on the output of the iVth syndes 
is no greater than -\/l\f/6* and thus that the rate_of accumulation 
of filtered waiting time jitter is no faster than \/N. Notice that to 
obtain this bound on the growth rate, it was not necessary to assume 
the stuffing ratios at each of the synchronizers identical nor the transfer 
functions of each of the desynchronizers identical. 

6.2 Experimental Data 

The power of the filtered waiting time jitter at the output of a chain 
of M12 syndes was measured for chain lengths N of 1, 2, 4, 8, and 16 
and six different stuffing ratios. The recorded powers are listed in 
Table IV and plotted in Figure 17. 

The upper bound + 

P«„v ^ N/6 

is plotted for comparison in each of the six graphs. The bound is typically 
about 15 dB above the data. No data are inconsistent with it. 
The other line plotted in each of the graphs is 

P».n = NP* nA , (16) 

where P* DA is the theoretical power in filtered waiting time jitter from 
just one syndes. The values used for P DA were taken from the computer 
program that produced the graph of Figure 9 and are listed in Table IV 
for convenience. This line was plotted because it seems to be a reasonable 
empirical approximation to the data. 

The fact that P u .i < Pu.\ for low stuffing ratios is rather surprising. 
We do not have a completely satisfactory explanation. For the stuffing 
ratios where this decrease in power occurs, the 1 "-lines and 1 Mines 



* It is of course being tacitly assumed here that the input to the first syndes is 
jitter free and that no other jitters (such as repeater jitter) are being introduced 
along the chain. A bound including the effects of other jitters could easily be obtained 
from (15) also. 

f There is about 0.15 dB peaking in the phase-locked loop of an M12 so we should 
strictly be using a bound based on (14) rather than (15). For this slight amount of 
peaking and the chain lengths considered here, the difference is negligible. 
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Table IV — Accumulation Data 



Stuffing 
Ratio 


Jitter Powers. 


in dB (0 dB = 1 Slot Squared) 




Theoretical 


Experimental 


p 


N = 1 


N = 1 


N = 2 


N = 4 


N = 8 


N = 16 


0.018 
0.092 
0.167 
0.241 
0.316 
0.391 


-11.8 
-17.3 
-21.5 
-22.9 
-22.6 
-24.5 


-10.9 
-14.8 
-20.7 
-21.8 
-20.5 
-20.7 


-13.1 
-20.7 
-18.4 
-19.6 
-15.2 
-19.8 


-10.0 
-17.1 
-19.0 
-16.8 
-14.9 
-15.6 


- 6.1 

- 9.5 
-16.8 
-12.8 
-11.7 
-12.5 


- 2.2 

- 5.7 
-10.3 

- 6.6 

- 9.1 

- 9.6 
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Fig. 17 — Experimental accumulation data as a function of log N. 
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of Q(/) are within the passband of //(/). It is felt that an explanation 
lies partly in the fact that Z t (j) may (and the approximation of Sub- 
section 5.1.3 suggests it will) throw much of the power in these lines 
out of the passband of #(/). If this idea is pursued and calculations 
made, however, it appears that while this phenomenon helps, it cannot 
be solely responsible. Multiple stuffing may possibly have played a 
role here also. 



VII. SUMMARY 

Expressions giving the spectrum of waiting time jitter both when 
there is and is not significant input jitter have been found. Using 
these theoretical expressions, graphs of the jitter power added by pulse 
stuffing versus the stuffing ratio have been drawn. These graphs indicate 
that when there is no input jitter there is much to be gained by intel- 
ligently choosing the stuffing ratio, but that with a typical amount 
of input jitter present, much of the advantage is lost. 

Bounds on the power in filtered waiting time jitter have been found, 
and it has been shown that when there is no peaking in the desyn- 
chronizers of a chain of syndes, the rms amplitude of filtered waiting 
time jitter accumulates at a rate no faster than the square root of the 
number of syndes in the chain. 
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APPENDIX A 

In this appendix the spectrum S s (f) of a (O is calculated. As noted 
previously, the procedure used will be to (i) write an equation describing 
waiting time jitter waveforms, (it) introduce initial condition random 
variables into this equation in such a way that a stationary ensemble 
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of waiting time jitter waveforms is defined, (in) compute the covariance 
of the waiting time jitter random process, and (iv) Fourier transform 
this covariance to obtain the power spectrum. 

Let [ • ] denote the greatest integer function, that is, let 



-1, -1 ^ x < 

r\ - 0, ^ X < 1. 

1, 1 ^ x < 2 



In Figure 18 a waiting time jitter waveform is drawn assuming stuffing 
opportunities occur at integer times and that just after stuffing at 
time t = 0, tp s (-) equals A — 1. With the aid of the greatest integer 
function an equation describing the waveform of Figure 18 can be 
written. It is 



4> s {t) = (A - 1) + pt - [p[t]]. 



(17) 



Obtaining this equation is the key step in finding the spectrum of 
waiting time jitter. The constant term. A — 1 is needed to make 5 (O) = 
A — 1. The second term pt generates the linearly increasing portion 
of the waveform <p 5 (0- Stuffs are made by [p[t]]. The reader should 
convince himself that this nesting of greatest integer functions puts 
stuffs in the proper locations. 

Equation (17) defines <t> s (-) at a stuffing time to be equal to its value 
just after a stuff has been made if one is to be made. In other words, 
it makes the function <f> s (-) continuous from the right. This convention 
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Fig. 18 — A waiting time jitter waveform starting at A — 1 at time t = 0. 
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is as good as any. For the purposes here it makes little difference how 
S (-) is defined at switching times. 

A waiting time jitter waveform drawn at random from a stationary 
ensemble of waiting time jitter waveforms (see Figure 19) will not in 
general (in fact almost surely will not) have stuffing opportunities 
coming at integer times. Define the random variable t as how long 
before time l = 0a stuffing opportunity last occurred. Over the ensemble 
t will be distributed uniformly on the interval [0, 1). 

Define the random variable f as by how much <t> s (-) exceeded A - 1 
at time t = -r + (see Figure 19). The random variable f must exist 
in the interval [0, 1). It cannot be less than because <f> s (-) never is 
less than (A - 1). It must be less than 1 because if 4> s (-) exceeds A 
just before the stuffing opportunity at time t = — r, a stuff will be 
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Fig. 19 — An ensemble of waiting time jitter waveforms. 
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Fig. 20 — Graphs of the functions u(x), v(x), and w(x). 



made. Over the stationary ensemble, f will be independent of t and 
distributed uniformly over [0, 1). 

A modified version of equation (17) that allows for nonzero f and t is 

<p s (t) - (A - 1) + r + Pit + t) - [f + p[t + r]]. (18) 

Define (see Figure 20a) 

u(x) = .v — [x] 

= x mod 1 . 

An equation equivalent to (18) is 

«/> 5 (0 = (A - 1) + u(T + p[l + r]) + pw(< + r). (19) 

The next step in the procedure is to find the covariance of the random 
process <f> s (-)- By definition the covariance C s (t) of <t> s (-) is given by 

C s (t) = E[(4> a (a + t) - ns)(4>s(s) - us)}, 
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where 

Us = E{Mt)} 

and E denotes expectation. 

The mean (t 8 is not difficult to evaluate. We have 

Ms = E {0,(0)} 

= E{(A - l)+1*(f + p[r]) + pu(r)\ 

= E\(A- l)+w(r)+ Mr)} 

= A - 1 + 1/2 + p/2. 

Define (see Figure 20b) 

v(x) = u(a) - 1/2. 

Then, for all i 

tpstt) ~ M.s = »(T + Pi* + r]) + pvit + r). (20) 

Returning to the evaluation of C s (t), we have 

C s (t) = E{(<p s (t) - m.s)(0,s(O) - i* s )) 

= E\v({ + p[* + r])i>(f + p[r])} + fl{e(f + p[t + r])pt>(r)} 

+ #{p!>(i + Mt + p[t])} + E[pv(t + t)p»(t)} 

= fl{w(f + Pit + r]Mf)} + E{v(S + p[* + r])p»(r)j 

+ E\pv(t + r)Kr)l + K\p<t + r)pw(r)}. 

The second and third terms in this expansion (cross covariances) are 
zero. Indeed, 

E[v(t + p[t + r])pv(r)} = P E\E\vU + p[< + t))v(t) I r}| 

= p E{v(t)E{v({ + p[t + r]) | r}} 

= P S{t;(r).0] 

= 

and 

£{pr(f + r)»(r)) = P E\E\v(t + r)Kf) | t\] 

= pE[v(t+ r)E{v(t) I r}} 

= P W* + 7-)^{«(r)}} 

= #{»(* + r)-0) 
= 0. 
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The fact that 

#{w(f + p[t + r]) | t) =0 

follows once it is observed that v(-) is periodic with period 1 arid that 
the integral of v(-) over any unit length interval is zero. Therefore, 

c,(t) - CsAt) + Cs. B (t), 

where 

Cs.aH) =E\vtt + p[*+ r]Mf)} 
and 

C,. B (t) = p 2 E[v(t+ t)v(t)\. 
The covariance C s .b(0 is the easier to evaluate. We have 
p- 2 C SlB (t) = E\v(l+ t)v(t)) 

= f v(t + t)v(t) dr 
Jo 

= f v{t + t)(t - 1/2) dr 
Jit 

= f ' (v(t) + t)(t - 1/2) dr 

+ f (v(t) + r - l)(r - 1/2) dr 

J 1 - ei ( ( ) 

= f (»(0 + t)(t - 1/2) dr - [ (r - 1/2) dr 

= f t(t - 1/2) dr - / (r - 1/2) dr 

J II Jl-u(t) 

= 1/12 - (1/2M/)(1 - u(t)). 
Define 

w(l) = 1/12 - (l/2)u(0(l - u(0). 
Then, 

Ca. B (0 = ptyfl. (21) 

A graph of «-(/) appears in Figure 20c. 

Turning to the calculation of C s ,a(0, we have 

c. s .,(0 -Jff(»(r + p[<+ r]MDI 

- J^l^Wr + p[< + rMf) |r}}. 
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The same manipulations that were used to find an expression for 
C a ,a(t) show 

E{v(t + pU + r]MD ! t) = w(p[t + t]). 
Therefore, 

C,,,(0 «£{w(p[* + 7-1)} 

= / w(p[t + r] dr 

= /" " w(p[t])dT+ f W(p[l+ l])dr 

Jo •'l-u(C) 



= (1 - u(t))w( P [t\) + u(t)w( P [t + 1]). 
A typical graph of C s , A (t) is shown in Figure 21. 



(22) 




Fig. 21 — The eovariance C.s\/i(0- 

It is convenient before calculating S s (f) to mention an equivalent 
formula for C s , A (t). Define 

1 - Ml* M I ^ 1. 
0, | * I > 1 

Then (see Figure 21), 

CbM = Z w(pn)A(t - n) 



A(t) = 



= A(t)*(w(pt) -rep 8(t)) 
where * denotes convolution and for any function X(j) 



(23) 
(24) 



repZ(/) = £ X(j-k). 
The Fourier transform of C s (t) is the spectrum of <p s (t). As a matter 
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of definition 

S s (f) = S[C 8 (t)} 

= 9[CmM] + SiCsAQ] 

= S s , A (f) + S S . B (J). 

The transform S s .n(f) is the easier to take. The function w(-) is 
periodic with period 1 and has zero mean. Thus, it has the Fourier 
series expansion 



It! 



(x) = 2 cj 2 ™ 



n» — oo 
n*0 



where for n = ±1, ±2, • • • 

c„ = / w(a;)e~ ,2irM d.r 

Jo 

= [ (1/12 - (l/2)u(x)(l - u(x)))e- i2rnx dx 
Jo 

= -(1/2) f u(x)(l - u{x))e~ i2rnz dx 
Jo 

= -(1/2) f x(l- 

Jo 

= {2imY 2 . 
Therefore, 

S S . B (J) = 5[C a . a (()} 



/7m 



-'.Ki)' 1 ^'"' 

= P 2 Z (^JW - ») + «(f + n)). 



(25) 
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Ss.aU) is only slightly more difficult to compute. Using equation 
(24) for C s , A (t), we have 

S a , A (f) = 5{C 8lA (t)} 

= 5{A(t)*(w(pt) -rep 8(t))} 

= 5[A(t)\-${w(pt) -rep 8(t)\ 

= ff{A(0}-(ff{w(pO}*ff{rep«(0}). 



The evaluations 



and 



_ , . . .. , sin x/ a . 2 , 
S[MW = -TZflT = smc / 



SF {rep 5(0} = rep 8(f) 



can be found in most Fourier transform tables. Proceeding as in the 
evaluation of S SiB (f), we have 

l n*0 > 

n*0 

= §(i) 2(i(/ - n ' ,) + w + ' , '' )) - 

Assembling these three evaluations, we obtain 

Ss, A (f) = sinc 2 /-Q(/), (26) 

where 

Q(f) = rep 8(f)* g (j^JW - *p) + (/ + »p)) 

= E (g^JW *tf - n p) + re P 5 (/ + "/>))• ( 2? ) 

APPENDIX B 

In this appendix the arguments of the previous appendix are gen- 
eralized to allow for input jitter. Rather than trying to write an equation 
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for 0s(O directly, it is easier to first write an equation for <pspc(t). 
With the jitter sign convention of Section V, the instantaneous 
frequency of the signal at the input to the synchronizer is / + 0/(0 
and the instantaneous stuffing ratio is p — 0/(0- As in Appendix A, 
let A denote the stuffing threshold, let t be a random variable denoting 
how long before time / = a stuffing opportunity occurs, and let f 
be a random variable denoting how far above A — 1, <t>spc(,') was 
just after the stuffing opportunity at time t = — r. The random variables 
f and t will again be independent and distributed uniformly on [0, 1). 
If there were no stuffing, 0s/-c(O would evolve as 

4>apc(t) = (A - 1) + f + / (p - 0,00) da. (28) 

In actuality there is stuffing, and at every stuffing opportunity a 
pulse is stuffed if 0s,-c(O is greater than A at that time. An equation 
describing 0s PC (O when stuffing is taking place is* 



4> S pS) = (A - 1) + f + f (p - 0,00) 
- [f + / (P - 0,00) ds] 



ds 



= (A - 1) + f + p(t + t) - 0,(0 + 0,(~r) 

- [r + P[t + r] - <p,([t + t]- t)+ 0,(-t)]. (29) 

A more convenient form for (29) can be obtained by again introducing 
the function u(-). It is 

*SPc(t) = (A - 1) + U(t + p[t + r] + 0,(-r) - 0,([/ + r] - r)) 

+ pu(t + t) + 0,([/ + t}- t) - 0,(0. 
Since 

0.,(O = 4>spc(t) +0,(0, 



* Actually, this equation stufTs more than one time slot if the waveform <pspc(^) 
immediately before a stuffing; opportunity has not only increased to a value above A, 
but is also above A + 1 (see Figure 22). A related problem arises if <f>.-ipc(l) immedi- 
ately before stuffing is below A — 1 (see Figure 23). The first situation will not 
occur if the sum of p and Ihe negative of the maximum possible negative change 
of 4>i(t) in one time slot is less than 1. The second situation will not occur if p minus 
the maximum possible positive change of #,(/) in one time slot is greater than zero. In 
other words, the model will not deviate from actuality if p is not close to or 1 and 
<t>i(l) is a low-frequency jitter. These conditions are reasonable and will be assumed 
to be met. 
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we have 

4> a (t) = (A - 1) + W(r + PU + r] - Mti + r] - r) + cf>,(-r)) 
+ pu(t+ r) +fa([t+ t] - t). 

The mean /x s of S (O is not difficult to compute. The integral of 
u(x) over any unit length interval is 1/2, and therefore 

E{ P u(t + t)\ = p/2. 

Similarly, 
Elutt 4- P [< + r] - *,([< + t]- t)+ *,(-t))} 

= ^!^|u(r+ p[«+ r] -*/([«+ r] - r) + 0,(-r)) | t, *,(■)} } 
-Jf|l/2| 

= 1/2. 



^iTtn^TT^^m^ 



*=S A 

^Ia-iI- — -- !l-^-K 

J L 



J I L 



J L 



i. i L 




Vwh 



1 I I L 




Fig. 22 — Deviation of the model from actuality when p = 1 . 
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Fig. 23 — Deviation of the model from actuality when p = 1 + . 



Therefore, 



and 



t, = 1/2 + p/2 



*«(0 - M.S = W(r + P[t + T] - *l([< + T] - T) + */(-T» 

+ pv(< + t) + (*,([« + r] - t) - M/) (30) 

where 

fJLj = A'!0/(O)- 

The three terms in this equation for <f> s (l) — m.s are easily shown 
to be uncorrelated. Thus, 
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Cs(t) = E{(4> s (t) - |U S )(0 S (O) ~ Ms)} 

= E[v(? + p[l + r] - MU + t\- r) + f(-r)Mf)J 
+ J5?M* + r)p»(r)]} 
+ 5{(*j([Z + t]- t) - M/)(0/(-r) - Mr)} 

+ C,.«(0 

+ c,.. t .(0. 

From Appendix A we have 

C s . B (t) =E\p\(t+ t)v(t)} 

= P 2 w(t) (31) 

where 

w(t) = 1/12 - (l/2)u(0(l - u(0). 

The evaluation of C s , A (t) is similar to the evaluation of the cor- 
responding term in Appendix A. We have 

Cs.aH) 
= E\v({ + p[t + r] - */([/ + r] - r) + fr(-r)Mf)} 
= £|£{y(r+ p[t+ r]-0/([/ + r] - r) +*/(-r))o(f) | r, 0,(-)}} 
= E\w(p[t + r] - 0,([< + r] - r) + fr(-r))}. 
Since 0/(0 is assumed stationary, 

C.,,. 4 (0 = E\w( P [t + r] - 0,([/ + t]) + 0,(0))} 

- tf{tf{w(p[* + r] -*,([* + r]) +0/(0)) 10,(0}} 
= (1 - u(t))E{w(p[t] - MM) + */(0))} 
+ u(t)E{w(p\t + 1] - 0/([/ + 1]) + 0/(0))} 

= A(t)*(r(t) -rep 8(t)), (32) 

where 

r(0 = E\w(pt -0/(0 +0,(0))}. 

Turning to the final term, we have again using the stationarity of 
*z(0 
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Cs.a(t) = B[fa(lt + r]) - m/)(0/(O) - M/)l 

= B[E\to([t + r]) - M/)(0/(O) " Hi) I */(■)!) 
= (1 - u(0)B{(*i(W) ~M/)(0/(O) -M/)i 
+ ll(fl*i(*i([* + ID -/*/)(*/(<>) -Mi)) 
= A(0*(B{(*/(0 - M/)(0/(O) - M/))-rep 5(0) 
= A(t)*(C(t)-vepm. (33) 

The spectrum <S. S (/) of the random process 0. S (O is the Fourier trans- 
form of C a (t). Let S s . A (j), S 8 , B (f), and S s . c (f) denote the Fourier 
transforms of Cs, A (t), C s .n(0, and C s .c(t). From Appendix A 

>s*,.„(0 = s|C,.„(0) 

= E («£rVw - ») + *(/ + »))• (34) 

„_1 \-TT71/ 

The next easiest Fourier transform to take is that giving S s ,c(D- 
We have 

fl-.cO) = fflCa.cWl 

= !F{A(0*(C,(0-repS(0)} 

= fF { .1(0 } • (fF { C 7 (0 1 * 5 Srep 5(0 }) 

= sine 2 /•(£,(/)* rep 8(f)) 

= sine 2 /-rep £,(/). (35) 

The evaluation of S SA (f) is the most difficult. We have 

Ss.aU) = 9[C a .M) 

- : T!.4(/)*(K0 -rep 5(/.))) 

- SF|A(0}-ff{r(0-rcp 5(01 
-sine 2 1-Q(f), 

where 

flfl) - 5!K0-rep5(0| 

= 3? | r(0)* rep 5(f). 

The periodic function w(a;) has the Fourier series expansion (see 
Appendix A) 
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W 



<»-±(azf*~- 



Therefore, 

JF{r(<)} - Z (2^) ${E\exp {j2im( P t - 0,(0 + 0,(0)) 



=-00 



= n Z (^J 2 *(/ - An)*ff{A'{exp {-ftm(*i(0 - 0,(0))}}} 

n»*0 

■ Jl GsJ *» - ** (36) 



where 
and 



Z n (f) = ff{z n (0} 



z„(l) = E{exp \-j2MM0 - 0/(0))}}. 
It is convenient to obtain another form for equation (36). Clearly, 
*-„(*) - E\exp {-JM-nKMt) ~ 0/(0))}} 
= tfjexp (-j2im(0 / (O) -0,(0)}}. 
Since 0,(-) is stationary, 

Z-M = ^(exp {-j27ro(0,(-0 - 0,(0))}} 

= *»(-*). 



Therefore, 



Z_„(/) 4 s\t-M\ 

= r 2 B (-06" ,ar/ ' dt 

= (" » li (0«""* ( " />l ^ 

J-oo 

= £„(-/). 
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Using this last relation in equation (36), we have 



Sir®) = E (^JV-tf ~ P") + Z-V + An)) 

= E (^W - pn) + Z»(-/ - pn)). 
Therefore, 

0(/) - rep S(0* E (^) >„(/ - P n) + Z.(-/ - pn)) 

= E (^Vep Z.tf - P n) + rep Z„(-/ - pn)). (37) 



In summary, 
S a (f) = sine 2 /•<)(/) 

+ E (2^) "(«(/ - ») + S(/ + n)) + sine 2 /-rep 5,(0 (38) 
where Q(f) is as above. 
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